Antiholons in one-dimensional t-J models. 
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Using a newly developed hybrid Monte Carlo algorithm for the nearest-neighbor (n.n.) t-J model, 
we show that antiholons identified in the supersymmetric inverse squared (IS) t-J model are clearly 
visible in the electron addition spectrum of the n.n. t-J model at J = 2t and also for J = 0.5t, a 
value of experimental relevance. 
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It is well established that electrons in one-dimensional 
(ID) metals generally lead to a Luttinger liquid 0,0, 
where charge-spin separation (CSS) takes place. The 
most direct evidence of CSS waspredicted for the spec- 
tral function of such systems 0, Q • While experimental 
evidence of CSS has accumulated in recent years H,IE0> 
an exact theoretical evaluation of the one-particle spec- 
tral function A{k, u>) could until now only be fully accom- 
plished for the Hubbard model at U = oo for arbitrary 
doping on the basis of the Bethe-Ansatz solution || |tj • 
However, recent progress was made for spectral proper- 
ties of the supersymmetric (SuSy) t-J model with 1/r 2 
interaction, where beyond the exact ground state (GS) 
[Toj . the thermodynamics the compact support of 
A(k,u)) 0| , the single-hole dynamics and the elec- 
tron addition spectrum |14j could be calculated analyt- 
ically. In addition to spinons and holons, the IS SuSy 
t-J model was shown to contain antiholons with charge 
Q = 2e, spin 5 = 0, and twice the mass of the holons, 
i.e. they are not merely charge conjugate to the holons. 

We present in this Letter one-particle spectral func- 
tions for the ID n.n. t-J model with finite doping, ob- 
tained by quantum Monte Carlo (QMC) simulations 
based on a newly developed hybrid algorithm. We show 
that the electron addition part of A(k,u>) presents a 
clear structure following the antiholon dispersion found 
in A(k,uj) of the IS SuSy t-J model for the same dop- 
ing. Our results show moreover, that also away from 
that point, a corresponding feature is present, strongly 
indicating that antiholons originally identified in the IS 
SuSy t-J model are generically present in the n.n. model. 

The t-J model reads: 
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Here. 



c\ are projected fermion operators 
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and cf and c i a are canonical creation and annihilation 
fermionic operators respectively with rii t(7 — c] G c i a . We 



consider two type of interactions: (i) n.n. type: = 
tSj,i+i, Jtj = J8j,i+i and (ii) IS SuSy type: i y = Jy/2 = 
t{-K I ' L) 2 1 sin 2 (7r(i — j)/L). L is the length of the system 
and the lattice constant is unity. For the n.n. model, 
it is convenient for the development of the algorithm to 
perform first a canonical transformation 
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where f\ , fi are canonical operators for spinless fermions, 
o ± = l/2(a x ±ia v ), and 7 ± = 1/2(1±ct 2 ), wither", a = 
x,y,z Pauli matrices. The Hamiltonian Q becomes: 

n = +tJ2 Pijflfj + 7 E - 1). (3) 
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where < i,j > means n.n. Pj, 
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and 



fjfj- The constraint against double 

occupancy becomes $^(1 — ff )// fi — 0, and commutes 
with the Hamiltonian. We consider now the following 
definition of an expectation value: 



(6) 



lim 



E n (^|Pe-T«Qe-f H ?|$ n ) 
E„(*n I V e-«« V | *„) 



(4) 



where | $f n ) =| s n ) ® | ^t), with {| s n )} a complete set 
of spin states and | ^t) a trial wavefunction for the spin- 
less fermions. V is a projector ensuring the constraint 
against double occupancy. Taking the limit O — > oo leads 
each state V \ #t) ® I s n ) to converge to the GS as long 
as the GS has a finite overlap with it. The multiplicity is 
corrected by the normalization factor. Introducing after 
slicing in imaginary time (typically time slices At = 0. 1/t 
are used) a complete set of spin states, and checkerboard- 
ing, the spin states are represented by world-lines, and for 
each configuration of them, fermions are evolved exactly 
since the Hamiltonian @ is bilinear in fermions. It can 
be easily shown that the total weight for a given config- 
uration of the world- lines is given by WnDf, where Wh 
is the weight of an antiferromagnetic Heisenber g m odel 
(AFHM), whereas Df is a fermionic determinant |l7l.ll8|. 
The updating of spin world-lines is performed using the 
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loop-algorithm [l9|. with the same complexity as for an 
AFHM, in contrast to a recently proposed pure loop- 
algorithm 20]. In fact, the autocorrelation time (r ~ 2 
for L = 30 and J/t = 2) for the internal energy is very 
similar to the one for the AFHM. Due to the mixed char- 
acter we denominate the whole hybrid-loop algorithm. 
Figure n shows a comparison of GS energies from QMC 
and exact diagonalization for various values of J at a 
density n = 0.9. The correct value is reached for values 
of the projection parameter 9 ~ 10/t — 20/t, demon- 
strating that the algorithm leads to the correct GS with 
high accuracy (statistical errors are smaller than the size 
of the symbols). Dynamical data arc obtained from the 
imaginary time Green's function and analytically contin- 
ued using the maximum entropy method |2l| . Further 
details will be presented elsewhere [T^ |. 
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FIG. 1: Ground-state energies vs. exact diagonalization re- 
sults for L — 20 with two holes as a function of the projection 
parameter O. 

We consider A(k,ui) both for electron-removal (ER) 
and electron-addition (EA) processes. Figure [21 a) shows 
A(k,uj) obtained from the QMC simulations for J = 2t 
(L = 40, 9 = 24/i), at a density n = 0.6. The Fermi en- 
ergy is taken as the zero of the energy scale. A splitting 
of the spectral weight into two branches can be read- 
ily seen on the EA side, in contradiction with what is 
expected for a single band. Figure [21 b) shows the pro- 
jection of A(k, w) on the {lo, k) plane, revealing the dis- 
persion of the main features in the spectrum, together 
with the compact support for the IS SuSy t-J model at 
the same density. Furthermore, the dispersions of spinon 
(s), holon (h), and antiholon (h) branches that deter- 
mine the compact support for EA processes in the IS 
SuSy t-J model are also shown. The dispersions for right 
(R) and left (L) going spinons and holons are given by 

£ S R(L)(?)A = <l( ±v s ~ ?)) and e hR(L)0?)A = Q(Q ± v c): 

respectively, where v® — w(l — n) and v® = it. The anti- 
holon dispersion is e^(q)/t = q (2v® — q) /2. The accessi- 
ble range of momenta is for and thR(L) i < q < k p 
( — kp <q<0), and for e s , < q < 2n~4k F ^M^M 
The compact support is obtained by assuming that the 
energy and momenta of the particle (EA) or hole (ER) 
are given by the addition of energy and momenta of s, 
h, and h with the dispersions above 01 • I n the ER part 



of the spectrum, only the corresponding part of the com- 
pact support and the dispersion of an antiholon branch 
along the support is shown for clarity, since in contrast to 
the EA processes, where only one spinon, one holon, and 
one antiholon are present, in the ER part three spinon 
and holon contributions [Q are possible. Therefore, a 
large number of features would appear, whose intensity 
is at the moment unknown, and hence, their importance 
is difficult to assess. A sharp feature is visible on the 
ER side that escapes from the compact support of the IS 
SuSy model. It is due to a holon branch, and as already 
discussed in the limit of a single hole |22|. the actual 
dispersion of the holon is needed, in order to describe 
this feature correctly. Also a deviation from the IS SuSy 
compact support is observed on the EA side at high en- 
ergies, where the differences in the models is expected to 
become noticeable. There are however, several features 
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FIG. 2: a) A(k,u) for J = It at a density n = 0.6. b) Pro- 
jection of intensities on the (w, k) plane. Solid lines: compact 
support of the IS SuSy t-J model. Red crosses: spinons, blue 
asterisks: holons, magenta diamonds: antiholons. See text 
for the dispersions. 

that are well described by the excitations of the IS model. 
The strongest feature on the EA side is followed closely 
by the spinon and holon branches between kp and 2kp, 
and for k > 2k f by a spinon at kp together with a dis- 
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persing antiholon. The analytic results for EA processes 
in the IS SuSy model [14j show that the largest portion 
of spectral weight is along this line. More striking is a 
second, weaker, but clearly visible branch that follows 
very closely the dispersion of an antiholon between kp 
and 2tt — 3&F- The analytic results of A(k,u>) for the 
IS SuSy model 0, |2^| predict a stepwise discontinuity 
at this edge and, in fact, the explicit evaluation of the 
weight shows for the present range of doping a higher 
value than in the interior of the support. Also the upper 
edge of the compact support on the ER part, that in the 
IS model corresponds to an antiholon, is well reproduced, 
with spectral weight down to very low energies around 
3fcp, as predicted by the IS SuSy model. Therefore, at 
the SuSy point, the clearest signal of CSS in A{k,u>) for 
the n.n. t-J model are present in the EA part of the spec- 
trum and through the comparison with the IS t-J model, 
it is clear that a sizeable part of the spectral weight goes 
to the antiholon excitation. Further analytical results 
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FIG. 3: a) A(k,w) for J = 0.5* at n = 0.6. b) Projection of 
intensities on the {u>, k) plane. Symbols coded as in Fig-HJ 

concerning the detailed structure of A(k,uj) on the EA 
part and for the one-holon and one-spinon contributions 
on the ER part will be published elsewhere |2^| . 

Since the exact solution of the IS model is restricted to 
the SuSy point, it is of much interest to see whether the 



features discussed above correspond to a generic behav- 
ior of the n.n. t-J model or whether it is better described, 
e.g. by the solution of the n.n. model at J = 8]. Figure 
a) shows A(k,uj) for n = 0.6 and J = 0.5* (L = 40, 
O = 56/*), i.e. very far away from the SuSy point and 
at a value of J/t of experimental relevance for cuprate 
compounds. A perspective was chosen, so that it is al- 
ready visible that as in the SuSy case, a structure splits 
off the main feature for k between 2kp and ir. Figure |3| 
b) shows the projection of A{k,u)) on the (w, k) plane. 
As a model for free spinons, holons, and antiholons, we 
use the same dispersions as for the IS SuSy model, but 
with e S R(L)(<z) = (^/2)?(±Wg — q), i.e. assuming that away 
from the SuSy point, only the energy scale of spinons is 
changed. The corresponding compact support, spinon, 
holon, and antiholon dispersions are encoded as in Fig. 
|3 In the present case, the compact support encloses 
rather well all the spectral weight. Moreover, on the ER 
part, the strongest feature is very accurately followed by 
a spinon, whereas a second structure is also closely fol- 
lowed by a holon. A more detailed view of these struc- 
tures is given below in Fig. They correspond to the 
generally expected signal in photoemission for CSS, that 
were also found in previous numerical studies of the Hub- 
bard model [24|]. However, as shown in Fig.0] the present 
algorithm seems to lead to results accurate enough, so 
that after application of maximum entropy, CSS is seen 
below Ep m. a, wider range in fc-space than previously. 
On the EA side, the feature with largest intensity is fol- 
lowed close to kF by a holon, a spinon, and an antiholon. 
However, further away from kp, the dispersion of the 
maximum is, up to k ~ 2kp, closer to an antiholon going 
from kp to 2ir — Zkp and beyond 2kp by a curve corre- 
sponding to a spinon at kp and a dispersing antiholon. 
Moreover, a second maximum develops beyond 2kp that 
follows the antiholon dispersing from kp to 2ir — 3kp, in a 
similar way as for J — 2t but with a smaller gap between 
both curves. In particular, the results from the simula- 
tions show appreciable weight between Zkp and the zone 
boundary, where only antiholons are present. 

A closer look to both features signaling CSS is given in 
Fig-El Figure0]a) shows A(k,uj) on the ER side and the 
location of the excitation energies for one spinon and one 
holon. Whereas the spinon dispersion follows the QMC 
data very closely, a deviation is seen for the holon for 
the farthest points from kp, as can be expected, since at 
higher energies, details of the dispersion matter in gen- 
eral. Yet, the agreement is good enough to enable an 
identification of the excitation content of the spectrum. 
The details of the splitted maxima for 2kp < k < ir on 
the EA side are shown in Fig. |Ub), where both an anti- 
holon dispersing from kp to 2ir — 3kp (closer to u> = 0) 
and an antiholon dispersing from 2kp to 2ir — 2kp on 
top of a spinon at kp are shown. Whereas the latter fol- 
lows the larger maximum, the former can be associated 
with the second maximum. As at the SuSy point, there 
seems to be almost no weight associated with the left 
propagating spinon and holon that give rise to the con- 
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tributions between 2kp and Zkp. This is consistent with 
the analytic results obtained for the IS SuSy model [l4| . 
Results for other values of doping (0.6 < n < 0.9) and J 
(0.5 < J/t < 3) not presented here, show the same qual- 
itative behavior, in particular the presence of a branch 
on the EA side below the main dispersing structure, that 
is closely followed by an antiholon branch under the as- 
sumption of free spinons, holons, and antiholons. 
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FIG. 4: Detailed view of A(k, u) at J = 0.5t and n = 0.6. a) 
ER side for < fc < fcf at J = 0.5i and n = 0.6. Vertical 
bars denote the positions for free spinon (closer to uj = 0) and 
holon excitations, b) EA side for 2&f < k < tt. Vertical bars 
denote here antiholon dispersions. 



In summary, one-particle spectra for electron-removal 
and addition were obtained using a new algorithm that 
delivers accurate dynamical data for the nearest-neighbor 
t-J model. A comparison with the compact support and 
excitation content of the 1/r 2 t-J model at the super- 
symmetric point J = 2t shows that a new manifestation 
of charge-spin separation in the n.n. model can be ob- 
served in the EA part of the spectrum, where in addition 
to spinons and holons, a branch following the antiholon 
dispersion is clearly visible. The same feature is still visi- 
ble at J = 0.5i, where assuming the same dispersions for 
the holon, and antiholon, as in the IS model but changing 
the scale of energy to J for the spinon, a fairly good de- 
scription of the spectrum can be given. Instead, serious 
deviations result by omiting the antiholon or setting its 
mass equal to that of the holon (i.e. assuming that it is 
the charge conjugated counterpart of the holon) |l8|]. The 
results above strongly indicate, that antiholons, that are 
not charge conjugate of holons, are generic excitations in 
the nearest neighbor t-J model. 
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